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Abstract. The dual condensate is a new QCD phase transition order parameter, which connnects confinement and chiral 
symmetry breaking as different mass hmits. We discuss the relation between the fermion spectrum at general boundary 
conditions and the dual condensate and show numerical results for the latter from unquenched SU (3) lattice configurations. 



INTRODUCTION 



DUAL CONDENSATE 



The QCD phase transition manifests itself in two phe- 
nomena, deconfinement and chiral symmetry restoration. 
The conventional order parameter for (de)confinement is 
the Polyakov loop, the straight loop in the time direction. 
After a suitable renormalization, it is related to the free 
energy of a static quark via (trP) ^ e^^^ (j3 = l/ksT). 
The Polyakov loop is small in the confined phase (and 
exactly vanishes in the quenched case because of center 
symmetry) and increases above the critical temperature. 

The chiral condensate, on the other hand, is an order 
parameter for chiral symmetry breaking in the massless 
limit, as it is not invariant under chiral transformations. 
In the chirally broken phase, the chiral condensate is 
finite, while it decays above the restoration temperature. 

Lattice simulations at physical quark masses have re- 
vealed the QCD phase transition to be a crossover with 
pseudo-critical temperatures of 157±4 MeV for the chi- 
ral susceptibility and 170 ±5 MeV for the Polyakov loop 
in [1] (see also [2]). 

The dual condensate [3] connects Polyakov loop and 
chiral condensate as two different mass limits, thus it 
also relates confinement and chiral symmetry breaking. 
It is therefore particularly interesting to see what one 
can learn from the dual condensate about the physi- 
cal crossover. We here improve previous results [4] and 
study the dual condensate on the Nf = 2 + 1 staggered 
dynamical configurations of [1]. 
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The physical boundary condition for a fermion field at fi- 
nite temperature is anti -periodic: \//(f + j3 = — \//(f ,x). 
With 'quark condensate' we refer to the expectation 
value E(m) = ^{Tr[{m + D_)^^]) with this boundary 
condition. 

We here also consider general boundary conditions [5] 



\l/it+P,x) =e'f\l/{t,x) 
giving the general quark condensate 



E(m,(p) = -(tr[(m + D^)-i]) 
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I ± iX,r. 
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where iX,p are the eigenvalues of the massless Dirac 
operator with these boundary condition (the physical 
boundary condition ^ = ;r is among them). 



FIGURE 1. Examples of closed loops on the lattice (with 
time running upwards). The red links get e'1' factors from the 
implementation of general boundary conditions, Eqn. (4). The 
green lines have winding number one. 

The dual condensate is defined as the first Fourier 
component of the general quark condensate with respect 
to the boundary phase (p: 
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FIGURE 2. Accumulated contributions of eigenvalues to tlie 
quark condensate Z()!j,;r) (top) and tire dual condensate Z(ra) 
(bottom), both in GeV^, at T = 172 MeV and m = 100 MeV. 
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FIGURE 3. Distribution of the lowest eigenvalues for the 
confined phase {T = 78 MeV, top) and the deconfined phase 
(T = 892 MeV, bottom). We compare histograms for (p = 
(red dashed) and (p = n (blue). 



The interpretation of this quantity is simplest in a lat- 
tice context. One can implement the boundary conditions 
(1) by multiplying a factor e'f to temporal links in one, 
say the last, time slice 

Uo{t = N,a) e'I'Uoit = N,a). (4) 

The general quark condensate E(m, (p) is gauge invari- 
ant and as such is composed of the contributions from 
all kinds of closed loops. These loops receive different 
powers of e'f factors, see Fig. 1. Now the dual conden- 
sate as the first Fourier component of E(ot, (p) (see (3)), 
picks out the contributions from all the loops with one 
e'f factor. These are loops winding once in the temperal 
direction, hence the dual condensate can be viewed as a 
'dressed Polyakov loop'. 

In a similar way dual observables can be constructed 
for arbitrary gauge invariant objects (cf. [6]). 

The conventional infinitely thin Polyakov loop is in- 
cluded in the set of loops that wind once and dominates 
in the limit of large probe mass m (which can be seen 
through an expansion in 1 /m). In this limit, however, 
more UV eigenvalues contribute to the sum in (3) [3]. 

We here consider all quantities unrenormalized (in [7] 
we demonstrated that the (quenched) dressed Polyakov 
loop has only a mild dependence on the lattice spacing). 

NUMERICAL RESULTS 

We use dynamical improved staggered fermion config- 
urations from [1] on lattices of size 8 x 24^, for tem- 
peratures ranging from 78 MeV to 890 MeV and lat- 



tice spacings from 0.282 fm to 0.028 fm. We compute 
500 to 1000 lowest eigenvalues of D for 16 or 8 different 
boundary conditions ^ G [0,2;r] with ARPACK. We cur- 
rently have completed the spectrum calculations for 20 
to 35 configurations at temperatures between 100 MeV 
and 200 MeV. 

The first problem we investigate is the convergence of 
the sums (2) and (3) when truncated to the number of 
available eigenvalues (in physical units). As the contri- 
bution of a ±/A pair to the condensates is 2m/X^ +m^, it 
is clear that this contribution decays for X ^ m and only 
the lowest part of the spectrum contributes. 

For the dual condensate there is an additional effect 
because it only probes the difference in the response of 
the spectra to changing boundary conditions. A strong 
response is manifest only in the IR spectrum [8, 6], as can 
be seen in Fig. 3, and for dual condensates thus only the 
IR contributes. Fig. 2 illustrates this effect: when using 
the available spectrum, the physical chiral condensate 
from (2) has not converged, while the dual condensate 
has due to the additional Fourier transform in (3). 

Fig. 4 shows the (unrenormalized) general quark con- 
densate as a function of the boundary angle (p at different 
temperatures. It is flat at low temperature and depends 
strongly on (p for high temperatures. Similar results were 
found also in non-lattice approaches [9]. 

Correspondingly, the dual condensate is small at low 
temperature and larger at high temperatures. It should 
serve as an order parameter for deconfinement, as the 
Polyakov loop does. In the quenched case this statement 
can be made exact because of the same behavior under 
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FIGURE 4. The general quark condensate (p) in GeV^ 
as a function of temperature and boundary angle with m = 
1 MeV (top). We remark that for large T and (p ^ n we expect 
further correlations until full convergence. The lower panels 
zoom into the confined phase (left, 7 = 78 MeV, m = 100 MeV) 
and the deconfined phase (right, T = 740 MeV, m = 10 MeV), 
respectively (each for a single configuration). 



center transformations [3]. Here both quantities have at 
least the same qualitative behavior. 

In Fig. 5 we show our results for the absolute value of 
the unrenormalized dual condensate as a function of tem- 
perature, compared to the conventional Polyakov loop. 
We also plot the negative logarithms of both divided by 
the inverse temperature (Fig. 6). For the Polyakov loop 
the latter has the interpretation of the free energy of an 
infinitely heavy quark. In analogy to that we might view 
the same quantity from the dressed Polyakov loops with 
mass parameter m as the free energy of a test quark with 
finite mass m. 

All of these quantities show an order parameter be- 
havior in the temperature range of 100 to 200 MeV. In 
the future we want to identify the critical temperatures 
(through inflection points and susceptibilities) and study 
their mass dependence. 

We thank Szabolcs Borsanyi for useful correspon- 
dence. F.B. and B.Z. are supported by DFG (BR 2872/4- 
2). 
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FIGURE 5. The dual condensate £(m) in GeV^ at m = 
60 MeV (top) and the Polyakov loop (bottom) as a function 
of temperature. 
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FIGURE 6. The 'free energy' — logZ/jS from the dual con- 
densate at m = 60 MeV (top) and from the Polyakov loop 
(— log |(tr -P) l/jS, bottom) vs. temperature. 
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